Exact solutions for a Dirac electron in an exponentially decaying magnetic field 
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We consider a Dirac electron in the presence of an exponentially decaying magnetic field. We 
obtain exact energy eigenvalues with a zero-energy state and the corresponding eigenfunctions. We 
also calculate the probability density and current distributions. 



PACS numbers: 03.65.Pm, 73.20.-r, 81.05.Uw, 71.70.Di 



The experimental realization of graphene, two- 
dimensional (2D) sheet of graphite p], and of the 
massless Dirac nature of its electron low-energy spectrum 
[3, 0] has given rise to a tremendous interest in this field 
(for recent reviews, see [f| @). The energy spectrum 
which goes linearly with the momentum and the specific 
density of states of the Dirac electrons 0] enabled the 
study of experimentally chiral tunneling and the Klein 
paradox in graphene [8(. This also leads to the anoma- 
lous Landau level spectrum in a uniform magnetic field, 
which gives rise to the half-integer quantum Hall effect 

The discovery of the half-integer quantum Hall effect 
and the zero-energy Landau level 0] stimulated lot of the- 
oretical research interest on Dirac electrons in uniform 
as well as non-uniform magnetic fields. For example, 
the Dirac- Weyl equation has been solved numerically for 
a single electron in a step-like magnetic field, magnetic 
barrier and polytropic magnetic field with B — B y' y , 
where 7 > [HJ EH, EH, E3j- We should also mention 
that the Schrodinger equation has been solved numeri- 
cally for spinless 2D electrons in a linearly varying mag- 
netic fie ld 14l| . step-like magnetic field, magnetic barrier 
US 0, 1 1 7t] - An analytical solution has been given for 
an electron in the presence of an exponentially decaying 
magnetic field [l8[ . 

In this work, we solve analytically for the Dirac elec- 
tron in graphene in the presence of an exponentially de- 
caying magnetic field B = Boe~ Xx . 

The Hamiltonian of the massless electrons in graphene 
near one of the Dirac K points is described by a two- 
component Dirac- Weyl equation 



H = vpo ■ p, 



(1) 



where vp w 10 6 m/sec is the Fermi velocity, a = 
{a x ,a y } is the 2x2 Pauli matrices and p = — ihV is 
the two-dimensional momentum operator. The pseudo- 
relativistic dispersion relation with Fermi velocity arises 
is due to the sublattice structure: the basis of the 
graphene honeycomb lattice contains two carbon atoms, 
giving rise to an isospin degree of freedom. In the pres- 
ence of an external magnetic field (B = V x A) perpen- 
dicular to the graphene plane, the Hamiltonian of the 
single Dirac electron (— e) is H = Vpcr ■ (p + eA). 

We consider an electron in the graphene sheet in 
the presence of a non- uniform magnetic field B(a;, y) = 



{0, 0, B$e Xx z} perpendicular to the x — y plane. In 
the Landau gauge, the corresponding vector potential is 
A(x,y) = {0,-{l/\)B o e- Xx y,0}. When A -> 0, the 
non-uniform magnetic field becomes a constant magnetic 
field. The time-independent Dirac- Weyl equation is 



v F a ■ (p + eA)*(a;, y) = E^(x, y). 



(2) 



Here, ^(x,y) 



{^+(2;, y), ^-(x, y)} T is the two- 



component wavefunction and T denotes the transpose 
of the column vector. Due to the translation invariance 
in the y-direction, the longitudinal momentum k y is a 
conserved quantity. One can parameterize the wavefunc- 
tion as ty(x,y) = e lkyy {ip+(x), , ip-(x)} T . From equation 
([2]) , one can get two coupled equations for ip + and ip- as 
given below: 



ifivp + {k y + ^Ay)]ip- (x) = Eijj + (x) (3) 



and 



ihvi 



dx 



(k y + -A y )]4, + (x)=E^{x). (4) 



Using two coupled equations. ([3]) and (|4]), we get 
Schrodinger-like decoupled equations for V>+ and ip- as 



o2 

t-(^) 2 ^2 + V +( X M'+ = ( 5 ) 



and 



d 2 



["(^) 2 ^2 + V-WW- = (6) 
Here, the effective potentials V± are given by 

V ± = (Kv F f[±~A y (x) + (k y + ^Ayf). (7) 

Here, the first term on the right hand side of the above 
equation is the Zeeman-like term due to the isospin de- 
gree of freedom in presence of the magnetic field. Define 
the magnetic length scale as l(x) = \Jh/eB{x). 

As seen in figure 1, the effective potentials V±(x) have 
the form of an asymmetric quantum wells formed by the 
exponentially decaying magnetic field. It is well known 
that such a well can have a bound state if the well is 



2 



sufficiently deep. Figure 1 will be discussed in more detail 
later on. 

Following the references. [H, [n|, we introduce two 
new dimensionless variables: 



1 eBne-^ 



(l(x)X) 2 h\ 2 



and 



\k y \_ eB e-^« 
SO = -r- 



A ?iA 2 («(a;o)A) i 



(8) 



(9) 



Here, 1(xq) = wh/eBjxo) is the magnetic length scale for 
a given value of x = xq which depends on the conserved 
wavevector k y through equation (8). Also, £ varies from 
to oo when x varies from oo to — oo. 

In the new variables, equations. (JSJ) and © reduce to 



1 d 



2£o Tl 



i]^± = o, 



(10) 



where /3 2 = £g — e 2 and e = E/{%vpX). The behavior for 
small and large £ suggests that the general solution can 
be written as ^±(£) ~ £,^s~^w±(£ l ). Inserting this ansatz 
in the previous equation, we get for w±(z = 2£) 



r d 2 , , d 

[z dV 2 + {l ~ z) Tz 



a±]w±{z) = 



(11) 



which has the form of a confluent hypergeometric equa- 
tion, where 7 = 2/3+1 and a± = — £o+/3+i±i. Two lin- 
early independent solutions can be chosen as F[a±,j; z) 
and U[a±,j;z], so the general solution can be written 
as w±(z) — AjF[a±, 7; z] + A 2 U[a±, 7; z]. Here, F and 
[/ are the first and second kind of confluent heypergeo- 
metric functions, respectively. However, U is not regular 
at the origin and has to be discarded. The requirement 
of normalizability implies that the solution is acceptable 
if a± is negative integer a± — —v^v — 0,1,2,.... This 
constraint produces the quantization of the energy: 



£4 



V(2£ ) 2 -(2£o-(2^ + l±l)) 2 . (12) 



The energy eigenvalues are then conveniently written as 



E n = hv F \y/(&) 2 - (£„ - nf 



(13) 



For E + , v = n — 1 and for E-, v = n. For n = and 
E + , v = — 1, but v can not be negative. For n = and 
£L, v = 0. Therefore, the n = state is not degener- 
ate. For n = 0, the solution u> + does not exist. We will 
consistently incorporate this fact by defining w n -% = 0. 
The corresponding wave function is given by 



w±(20=iFi[-(n- 



I 2/3 + 1.2& (14) 



where 1-F1 [— n, a, x] is the confluent hypergeometric func- 
tion. 



Equation (|13|1 can be re-written as 



E 2 l = h 2 uo 2 c (x )2n(l---) 1 



n 

Wo' 



(15) 



where the local cyclotron frequency is u c (xq) — 
{vf)/1(xq). Note that when A — * (£0 — *• °°), the 
energy eigenvalues reduce to the well known relativis- 
tic Landau level structure for uniform magnetic field: 
E n = (hvF/l(x ))V2n- 

When A is very large, the effect of inhomogeneous mag- 
netic field vanishes and the Dirac electron feels no effec- 
tive potential. It behaves like a free particle. In this limit, 
the energy spectrum becomes E n oc n which agrees with 
the known result for the free Dirac electron [201 ] . 

The complete normalized wavefunctions can be written 

as 



ik y y 



y/2L y l(x ) 



(26) 



and 



ik y y 



^2L y l(x ) 



-\f3X e -£ e- 



(16) 



(17) 



where L y is the length of the system along the y axis, 
Co = v2 and c„>o = 1. Also, X = (x — xq)/1(xq), A = 
XI (xq), and w n is given by 



r[2/3- 



1]A 



r[2/3+l]r[2/3]r[n+l] 



iiq[-n,2/3+l,2£ e 



-AXi 



(18) 

We have also checked that the wavefunctions (fTB"]) and 
(fT7|) reduce to that of the Dirac Landau level for constant 
magnetic field (i.e. a Hermite polynomial multiplied with 
a Gaussian factor) when A — > 0. The phase factor i 
in the lower component wavefunction (|17p is obtained 
from equation ([J), which is crucial for calculating the 
probability current density. 

The effective potential ([7]) can be re-written as 



V± = (hu c (xo)) 2 [±e~^ +&(1 - e~ AX ) 21 



(19) 



For weak-inhomogcneity (A << 1), the effective poten- 
tials are almost symmetric around X = point. For 
strong inhomogeneity (A >> 0), there is a strong asym- 
metry of the effective potentials. The effective poten- 
tials for negative k y does not have any global minimum. 
Therefore, the bound state does not exist for negative k y . 
The effective potentials together with the energy eigen- 
values for various values of A are shown in figure 1 . Both 
the effective potentials get saturated to (Tiu> c (xo)) 2 £,o at 
large X. The zero-energy state (n = 0) is always lies in- 
side the potential V- but outside the potential V+ . How- 
ever, all other discrete energy levels (n > 0) are lying 
inside both the potentials V±, which is expected from 
the solution of the Dirac- Weyl equation. The number of 
energy levels decreases as we increases A. For example, 
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FIG. 1: (Color online) Plots of the effective potential V±(x) 
(in units of (hu} c ) 2 ) vs X together with the energy eigenvalues 
E n (in units of %u> c ) for A = 0.5 and A = 0.01 (only zero- 
energy level is shown in this case). 



there are five energy levels including zero-energy state 
when A = 0.5. On the other hand, number of energy 
levels is quite large when A = 0.01. There are a finite 
number of discrete energy levels for a given asymmetric 
parameter A. The total number of discrete energy levels 
including zero-energy level can be calculated easily from 
the condition that E\ < V±(X — > oo) and it is given by 
N = Int[^o] + 1. Here, Int[b] means the integer that is 
just smaller than b. In figure 1, only the zero-energy level 
is shown for the A = 0.01 case. 

The probability density distribution of Dirac elec- 
trons in a given level n is p n (x) — vl^wf (x)^^ (x). 
The probability density distributions p n for different 
values of n and A are shown in figure 2. The ve- 
locity operator that follows from the Heisenberg equa- 
tion is given by v = vf<J. The probability cur- 



rent distribution is J n (x) = ev F ^ n ^ (x)a y ^ n ^ (x) = 

-ievpi^ix^ix) - 4rL n) *(z)#i n) (z))- The prob- 
ability current distributions for various values of n and 
A are shown in figure 3. The zero-energy state does not 
carry any current, irrespective of nature of the magnetic 
field. When A — > 0, p n and J n are symmetric around the 
point X — 0. When A >> 0, p n and J n are strongly 
asymmetry as is expected from the effective potentials. 

In summary, we have obtained exact energy eigenval- 
ues and the corresponding eigenfunctions for a Dirac elec- 
tron in th e presence of an exponentially decaying ma g- 
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FIG. 2: (Color online) Plots of the probability density distri- 
bution p n (x) (in units of \/{L v l{xo)) vs X for various values 
of n with A = 0.5 and A = 0.01. 



netic field. We have also provided the probability density 
and current distributions for each band. 
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FIG. 3: (Color online) Plots of the probability current density 
Jn(x) (in units of evF / (L y l(xo)) vs X for various values of n 
with fixed A = 0.5 and A = 0.01. 



